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Abstract 

We consider orthogonal polynomials {Pti,n(x)}^Lo on the real line 
with respect to a weight w(x) — e~ NV i x ) and in particular the asymp- 
totic behaviour of the coefficients a n ,N and b n ,N in the three term 
recurrence xn n , N (x) = n n+ i tN (x) + b niN n ntN (x) + a ntN iv n -i iN (x). For 
one-cut regular V we show, using the Deift-Zhou method of steepest 
descent for Riemann-Hilbert problems, that the diagonal recurrence 
coefficients a n ,n an d b n , n have asymptotic expansions as n — > oo in 
powers of 1/n 2 and powers of 1/n, respectively. 



1 Introduction 

We consider the asymptotic behavior of the recurrence coefficients a n ,N and 
b n N in the three-term recurrence relation 

X7T n ,N(x) = 1T n+1 ^ N (x) + b n)N 7T n ^ N (x) + a n>N TT n ^ liN (x) 
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a grant from the Ministry of Education and Science of Spain, project code MTM2005- 
08648-C02-01. 
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for orthogonal polynomials with respect to varying exponential weights. 
Here n nt N is the n-th degree monic orthogonal polynomial with respect to 
a varying weight 

w N (x) = e- NV & 

where V is real analytic on R with lim — Vff ^ = +00. Moreover, V is 

assumed to be one-cut regular, which means that the equilibrium measure 
dfiy = ipv(x)dx associated with V is supported on one interval [a, b] where 
it has the form 

t^v(x)dx = yjb — x)(x — ajh(x)x[a,b] i x ) dx (1-1) 

where h is real analytic, strictly positive on [a, b], and in addition the inequal- 
ity U3U is strict for x G R \ [a, b). See e.g. [H El dH [IT] for the definition 
of the equilibrium measure and for more information on the one-cut regular 
case. 

Under these assumptions Deift et al. [7j proved that a n ^ n and b Utri have 
asymptotic expansions in powers of 1/n. Their approach is based on the 
Deift-Zhou method of steepest descent applied to the Riemann-Hilbert pro- 
blem for orthogonal polynomials of Fokas, Its, and Kitaev [12]. This method 
was first introduced in [9] and further developed in (HJ [Jj, [8] and many papers 
since then. 

The asymptotic result on the recurrence coefficients was considerably 
refined by Bleher and Its [21 Theorem 5.2] who showed for polynomial V that 
there exists e > and real analytic functions f2k(s), 52fc(s), k = 0, 1, . . ., on 
[1 — e, 1 + e] such that the asymptotic expansions 

00 

a n ,N ~ /o (£) + £ iV" 2m / 2m (£) (1.2) 

m=l 

/n + l/2\ * t 9m fn + l/2\ , . 

hold uniformly as n, N — > cxd with 1 — e < n/A^ < 1 + e. These 1/N 2 
expansions are used in [2] to prove the 1/N 2 expansion of the free energy 
(a.k.a. logarithm of the partition function or Hankel determinant) of the 
associated random matrix ensemble in the one-cut regular case, see also 

mi- 

The proof of (jl.2p and (jl.3p in [2] is based on the Deift et al. result 
referred to above, in combination with so-called string equations. It is of 
some interest to find a proof that is based on the Riemann-Hilbert steepest 
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descent analysis only. Here we do this for the diagonal case n = N, and we 
obtain the following. 

Theorem 1.1. Let V be real analytic and one-cut regular. Then there exist 
constants and (3 m , m = 1, 2, . . . (depending on V) such that a„ jn and 
b ntn have the following asymptotic expansions as n —* oo: 

(ft - a) 2 . v> «2m , b + a A /3 m 

On,n + 2^ V« ^~ + ( L4 ) 

m=l m=l 

where a and b are the endpoints of the support of ifty. The first coefficient 
(3\ in the expansion for b n<n is given explicitly by 

01 = 2vr(6 - a) (h(bj ~ %)) ^ 

where /i is the function appearing in the expression (jl.ip for the equilibrium 
measure ipy associated with V. 

In our proof of Theorem 11.11 we follow the main lines of the steepest 
descent analysis of [7] . We will deduce that the odd powers in the expansion 
of 

a n,n vanish from the structure of the local Airy parametrices around the 
endpoints. The expression (|1.5p for f3± is new, although it is likely that it 
can be deduced from the approach of [2] as well. The explicit formula (jl.5p 
shows that f3\ = if and only if h(a) = h(b). It is very easy to construct 
examples of one-cut regular V such that h(a) / h(b) and so f3\ ^ 0. We 
have thus corrected an error in a paper of Albeverio, Pastur, and Shcherbina 
[U Theorem 1, formula (1.34)] who claim that (3\ = always in the one-cut 
regular case. 

Example 1.2. We may explicitly check Theorem 11.11 using Jacobi poly- 
nomials with varying parameters a = AN, (3 = BN, A,B > 0. These 
polynomials are orthogonal with weight (1 — x) AN (1 + x) BN on [—1, 1]. The 
equilibrium measure takes the form (jl.ip with 



B 2 - A 2 ± V(l + A + B){l+A)(l + B) 
a,b= — : — t, (l-o 

(2 + A + B) 2 v ' 

and 

2 + A + B 
27r(l — x z ) 

see [H1[T5]. We are in the one-cut regular case, but for weights restricted 
to [—1, 1]. An analysis of the proof of Theorem ll.il however, will show that 
the results (ll.4|) - (ll.5p remain valid in this case as well. 
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From the explicit form of the recurrence coefficients for Jacobi polyno- 
mials, see e.g. [H [T5]. 

4(1 + A + S)(1 + A)(1 + S) 



((2 + A + B) 2 -^){2 + A + B) 2 

B 2 — A 2 

{2 + A + B){2 + A + B + lY 



it is easy to see that (|1.4p holds. Using (|1.6|) - (|1.7|) we can also ascertain the 
validity of (fT3|) . 



2 The Riemann-Hilbert Problem 

The Riemann-Hilbert problem for orthogonal polynomials was introduced 
by Fokas, Its, and Kitaev [12]. It asks for a 2 x 2 matrix valued function 
Y(z) satisfying 



Y(z) is analytic in C \ 
Y+(x) = Y-{x) ( „ , ) for x G 
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(2.1) 



y ( z) = (/ + 0(i)) (; n z i 



The unique solution of (|2.ip is (see e.g. [5]) 



Y(z) 



i 



as z — > oo. 



\ 



-2iriKn-i t NPn--L,N(z) -Kn-1,N J 



2mn n jv J t— 2 
' K 

Pn-l,jvft) 



t-Z 



(2.2) 



/ 



where p n ,N{x) = K n ^7r n ^(x) is the nth degree orthonormal polynomial. 
The recurrence coefficients are expressed as follows in terms of the solution 
of the Riemann-Hilbert problem (|2.ip . see [5| IIP]. 



Proposition 2.1. Let 

Y(Z): 

Then 



1 1 f I 

7 + - Y x + - Y 2 + - 

z Z L \ z° 



z n 
z~ n 



a n ,N = (Y x ) 12 (Yi) 
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(2.3) 



(2.4) 
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and 



(Yi) 



(2.5) 



(Yi) 
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For the remainder of this paper we will take N = n. We closely follow 
[5j [7j in applying the Deift-Zhou method of steepest descent for Riemann- 
Hilbert problems to (|2.1j) . 

3 The Deift-Zhou method of steepest descent 

The goal of the Deift-Zhou method of steepest descent for Riemann-Hilbert 
problems is to change the original problem into a problem for which the 
asymptotics for z — ► oo are normalised and for which all matrices, jump 
matrices and solutions alike, are asymptotically close to the identity matrix 
for large n which can be solved iteratively. The specific details and steps 
needed to achieve this goal shall be explained below. 

3.1 The First Step: Transformation Y i— ► T 

The key aspect of the first step of the analysis is the equilibrium measure \xy 
corresponding to V. This equilibrium measure [iy is the unique probability 
measure that satisfies for some I, 



2 / log \x — y\ dfiy(y) + V{x) = I, for all x G supp pLy (3-2) 



For the one-cut regular case that we are considering we have that the support 
is one interval [a,b] and dpLy{x) = ij)y(x)dx as in (jl.ip . In addition the 
inequality (|3.ip is strict for x G M \ [a, b]. 




for all x G E 



(3.1) 



Define 




(3.3) 



and 




(3.4) 



(3.5) 
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If we now put 



T(z) = e " n ( / / 2 ) fT3 y(z)e _n9 ^) f73 e n ^ //2 - )<T3 , 



(3.6) 



where 03 = ( j-j -1 ) ^ s the third Pauli matrix, then T satisfies the Riemann- 
Hilbert problem 



T(z) is analytic in C \ R, 
T + (x) = T^{x)Jt{x) for x € 1 
k T(z) = 1 + 0(1) as z^oo, 



(3.7) 



where 



^ e -2n<p(x) 

1 

e 2n<t>+(x) ^ 

g 2n0_(x) 
^ e -2n</)(x) 

1 



for x < a, 

for re G (a, 6), 
for x > b. 



(3. 



Since the inequality in (|3.1j) is strict for x < o and x > b we have that 
0(x) > for x < a and <^(x) > for x > b. Thus the jump matrices for T 
on (— oo, a) and (6, oo) tend to the identity matrix as n — > oo. 

3.2 The Second Step: Transformation T h-> 5 

The second transformation is the so-called opening of the lens and it is based 
on the factorisation 



,2n<j) + (x) 





1 

J2n$-(x) 



1 0W0 1W 1 

e 2n<t>_(x) 2/1—1 0/ l e 2n< M a; ) 1 



(3.9) 



of the jump matrix Jt on (a,b). The factorisation (|3.9p allows us to split 
the jump on (a, b) as shown in Figure [TJ 

We use Si and £2 to denote the upper and lower lips of the lens, respec- 
tively. We define S as follows: 

• For z outside the lens, we put S = T. 

• For z within the region enclosed by £1 and (a, b), 

1 (T 



S = T 



2nct> 1 I • 



(3.10) 
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For z within the region enclosed by Si and (a, 6), 



S = T 



1 

e 2n<P 1 



Then S satisfies the following Riemann-Hilbert problem: 
' S(z) is analytic in C \ (R U S x U S 2 ) 
< S+{z) = S-{z)J s {z) for z G R U S x U S 2 
k = / + Q) for z -> oo 



where 



e 2n^(z) J) for z G Si U S 2 , 
for z G (a, 6), 
for z < a, 
for z > b, 



(3-11) 



(3.12) 



(3.13) 



1 
-1 

1 g-2n<Mz) 

1 

1 g-2n</>(2) 

1 



We may (and do) assume that the lips of the lens are in the region where 
Re (f> < 0, so that the jump matrices for S on Si and S 2 tend to the identity 
matrix as n — > oo. 
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3.3 The Third Step: Parametrix Away From Endpoints 



The parametrix away from the branch points is a 'global solution' N(z) 
satisfying the Riemann-Hilbert problem 

N (z) is analytic in C \ [a, b] 

(3.14) 



N+(x) = N-(x) l_ 1 JforxG(a,6) 
k iV(z) = 1 + 0(1) for z^oo 
which has solution (see [5]) 



ffOQ-/?- 1 ^) /3( 2 )+ 2 /3~ 1 (z) 
2i 2 



(3.15) 



where /3(z) = fff 



-6 \ 4 



3.4 The Fourth Step: Parametrices Near Endpoints 

Having constructed the 'global solution', the next step is finding 'local solu- 
tions' close to the endpoints a and b. Near b, the local situation is described 
as in the left picture of Figure [2] with jump matrix 

f < 1 0" 

e 2n4>(z) i 

' 1 N 
-1 o y 

e -2n</>(«) N 
,0 1 



Jp(z) = Jsfc) = < 



on Si n U and S 2 n U 
on (a, 6) n U 

on (6,oo) n C7 



where U is a (small) disk around b. 

We therefore want to find a matrix function P, that solves 

P(z) is analytic on U \ (Si US 2 UR) 
P+(z) = pJ(z)J P (z) on (Si US 2 UR)n[/ 
P(z) = iV(z) (l + (^)) as n — > oo uniformly for z E dU 

Then p(z)e n ^ <yZ ^ (7 ' i should have constant jumps on (Si U S2 UR) n U, namely 

'l N 



P{z)e 



n4>(z)cr3 



P{z)e 



n4>{z)crz 



X < 



1 1 

1 
-1 

1 1 

1 



for z e (Si u s 2 ) n u 

for z € (a, b) D U 
for z G (6, 00) fl U 



s 




Figure 2: Mapping of neighbourhood of b onto a neighbourhood of f(b) = 



Shrinking U if necessary, we have that 
C = f(z) = (\<t>{ 



2/3 



defines a conformal map from U to a convex neighborhood of £ = 0. We may 
and do assume that the lips of the lens are taken so that Si n U is mapped 
into argC = 2ir/3, and £2 H U is mapped into arg£ = 27r/3, see Figure [21 
Denoting the sectors in the C-plane by I, II, III, IV as in Figured! and using 
the usual Airy function Ai(£), we construct the Airy model solution by 



*(C) 



/Ai(C) uAi(uC) 
\M'(0 W 2 AiVC) 
/Ai(C) -^ 2 Ai(^ 2 C) 
VAi'(C) -wAiVO 
-wAi(wC) -w 2 Ai(a; 2 C) 
-u 2 Ai'(uj() -toAi'{uj 2 C) 
-uj 2 A\(uj 2 C) uAi(u() 

-L0Ai'(LU 2 () LU 2 Ai'{cu() 



for £ in sector IV 
for £ in sector I 
for £ in sector II 
for £ in sector III 



which has the jump matrices in the £-plane indicated in the right side of 
Figure El 

Then for any analytic prefactor E n (z) we have that 



P(z) = E n (z)<5>(Jf(z))e n ^ 
has the required jump matrices Jp. If we choose 



n 



2/3 



/(*) 



0-3/4 



then the matching condition P(z) = N(z)(I+0(l/n)) as n 
is satisfied as well, see e.g. [31 [51 E] for further detail. 



(3.16) 

(3.17) 

00 for z £ dU, 
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A similar construction yields a parametrix P in a small disc U around a. 
One can see that P can be obtained by taking P and interchanging a and b 
and conjugating with 03. 

3.5 The Fifth Step: Transformation S \-> R 

Using the parametrices N, P, and P, we define the third transformation 
S 1 ^ R as follows 



5(z)7V(z)- 1 for z G C\ (UUU) 
slz)P(z)- 1 for z £ U 
s\z)p\z)- 1 for z £ U 



(3.18) 



Then i? has no jump on [a, b] \ (U U U), as the jumps of S and iV _1 cancel 
out. In U and C/ the jumps of S cancel out with the jumps of P and P, 
leaving only jumps for R on the contour Y,r shown in Figure El 




Figure 3: Contour S/j for the Riemann-Hilbert problem for R 
The Riemann-Hilbert problem for R is 

R(z) is analytic on C \ 
R+(z) =R-(z)J R {z) forzGSR 
R(z)=I + (l) for z^ 00 



where 

r N(z) J s (z)N(z)- 1 for z G S fi \ (at/ U 

J fl (z) = J P(z)iV(z)- 1 for zGdU 

{ P(z)iV(4 _1 for 2 G <9C7 

The jump matrices Jr(z) = N(z)Js(z)N(z)~ 1 tend to the identity ma- 
trix at an exponential rate as n — > 00. The jump matrices on dU and <9t/ 
tend to the identity matrix but at a slower rate of 1/n as n — > 00. The 



10 



precise form is obtained from the asymptotic expansion of the Airy function 
as z — ► oo, — 7r < argz < ir, (see [T3l) 



o 3 
- 2 r2 CO 



Ai(js) r > - , , r — ftt^^t- 3.19 

1 ' 2^fcZ 9*(2fc)!r(J) z i* 1 J 

and the corresponding asymptotic expansion for Ai (2). Using these facts 
in the parametrix P we find an asymptotic expansion for the jump of R on 
dU 

J fi (z) = P^iV^)- 1 ~ / + £ ^A fc (z) (3.20) 

where 



jfc=l 



a ^ = J_ / r(3fc + i) r(3fe-|) \ 1 

fel J 9*(2&)! 4-9*- 1 (2(*-l))V (§0(*))* 

1 r(3fc-|) 1 



V5F9*-i(2(fc-l))!(3 0(z) ) 



70"2 for even (3-21) 



and 



a , x Z 3 ^) 2 1 / T (3A; + i) r(3fc-f) V 



)) fe 20F\ v 9 fc (2/c)! 2 • 9 fe - 1 (2(A; - 1))! 

03 — i(7i) for k odd (3.22) 



/3(z)~ 2 1 r(3fc + i 



where 



(§0(2))* 2v/5F 9*(2*)! 



(v i? o')' »»=(J ,:! ' 2:il 













0"2 = f 



are the Pauli matrices. 
A similar expansion 



J fl (z) = P^)iV(z)- 1 ~ / + ^ — A fc (z) (3.24) 



fc=i 



holds for the jump matrix on dU. 

As a result we find by the methods of [7J, see also [14|, Lemma 8.3], 
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Lemma 3.1. There exist matrix valued functions Rk(z) with the property 
that for every I £ N, there exist constants C > and r > such that for 
every z with \z\ > r, 

l 

k=l 

So we write 

R(z) ^I + J2-^Rk{z) (3.26) 

k=l 

From (|3.26p . (|3.20p and H3.24H and the Riemann-Hilbert problem for R, we 
find an additive Riemann-Hilbert problem for R k (z), 

R k (z) is analytic on C \ (dU U dU) 

k— 1 

R k+ (z) = R k ^{z) + £ /2,_(z)A fc _,(z) for z G dU 

k-°t ( 3 - 27 ) 
i? fc+ (z) = R k _(z) + E i?,_(z)A fc _ { (z) for z e 8U 
l=o 

Rk(z) = O (-) as z — ► oo 

where Rq(z) = I. These Riemann-Hilbert problems can be successively 
solved using the Sokhotskii-Plemelj formula, or using a technique based on 
Laurent series expansion as in [14] . 



< 



C 



\z\n 



i+i 



(3.25) 



4 Proof of Theorem 11.11 

For the proof of (|1.4[) we do not need to compute the explicit forms of the 
R k s. However, we need to know that they have the following structure. 
Recall that the Pauli matrices are given in (|3.23j) . 

Lemma 4.1. For k odd, R k (z) is a linear combination of o\ and 03 and for 
k even, R k (z) is a linear combination of / and 02- 

Proof. For k = 1, we know because of (I3.27P that R\ + = + Aj on dU 
and = + Ai on dU. As Ai, Ai G span {<ri, 0-3} on account of 
()3.22p . R\(z) must be a linear combination of o~\ and cr 3 as well. 

Let k > 1 and once more observe (|3.27p . If is odd, then again by 
(|3.22p A k , A k £ span {<7i, 0-3} and using induction on k, for every I < k, 
Ri^(z)A k ^i(z) and Ri-(z)A k ^i(z) are products of a linear combination of 
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o\ and (T3 and a linear combination of / and 02 (see also (|3.21|) f)3.22[) ). 
which results in a linear combination of o\ and 03. Thus all terms in the 
(additive) jump for on dU and on dU are in the span of a\ and 03, and 
it follows that Rk £ span {01, 03} if k is odd. 

If is even, then by induction, where we use again ( (|3.21|) - (|3.22|) . we 
have that Ri-(z)Ak-i(z) and Ri-{z)Ak-i{z) are either products of two li- 
near combinations of I and 02 (in case / is even), or products of two li- 
near combinations of o\ and 03 (in case I is odd). In both cases we find 
that Ri-(z)Ak-i(z) and Ri-(z)Ak-i(z) are linear combinations of / and 02, 
which implies that R^ £ span {1,02} if k is even. □ 

Now we can finally prove our main result. 

Proof of Theorem We start from the expressions (12, 4D and (12. 5p for a n ^ n 
and b n<n in terms of the solution of the Riemann-Hilbert problem for Y . 
Following the transformations Y 1— > T 1— ► S, we find that 

a n ,n = (5i) 12 (5i) 21 (4.1) 

and 

= " (^l)22 (4-2) 

where Si and S2 are the terms in the expansion of S(z) as z — » 00, 

= J + + + O (J^ 

To obtain (14.2)) we use that g(z) = logz + (D(l/z), see also [TO] , 

By (|3.18p . we know that S(z) = R(z)N(z) for \z\ large enough, so we 
need the first terms in the expansions of N(z) and R(z) as z — > 00. From 
(|3.15p we have 

f3(z) + f3(z)-\ , 0(z)-0(z)-\ 
N{z) = 1 H cr 2 

r (b-a) 1 /(^-a) 2 r ^-a 2 \ 1 -fl\ , A n , 

= / - l x Jct2 ; + S^ 1 — r-* )? +o ta0 (43) 
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and from Lemma |4,1I 



R(z)=I+-( ^ {Rml<n°l +Rmla z ^2) + ^ {Rmlll + Rmia^) J 

\m odd m even / 

H j^f ^ (Rm2a 1 °'l + Rm2a 3 0'3) + ^ ~ ^ {Rm2ll + Rm2a 2 a 2) ] 

\m odd m even / 

+ (±) (4.4) 



where the constants R m ji, R m j <Tk , for m G N, j = 1,2, and k = 1,2,3 are 
such that Rmjil + Yl\=iRmja k &k is the coefficient of z - - 7 in the Laurent 
expansion of R m {z) around 2 = oo. 
Therefore, by (g3]) and (|Q|> . 

=R{z)N{z) ~ / + 1 ( - ^"^ CTS + V + ^mlos^) 

\ m odd 

+ ^ — ^ {Rmlll + Rmla 2 <T 2) J 
m even / 



1 /(&-a) 2 & 2 -a 2 iff b - a \ 

* I —32— 7 " — §— ^ + 2^ ^^m2 CT1 + i?ml CT3 J 0"! 



, .b — a 

+ I Rm2cr 3 — I — I — 



+ ( Rm2a 2 ~ ^-T^Rmll) <?2 ]) + O ( \ ) (4.5) 



which implies that 

(Si)i2 ~ h —r L i + X) i^-i-i " * X ^ Rm ^ ( 4 - 6 ) 



m odd m even 



and 



m odd m even 



Inserting (|4.6p and (|4.7p into (|4.ip then finally gives 

(6 - a) 2 ^ a2r 



- — — — i- 



16 ' ^ n 2m 

m=l 
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for certain constants ai m . 

Similar to (|4.6p and (|4.7p we have that (52)12 and (£1)22 have asymptotic 
expansions in powers of 1/n. Prom the expansion (|4.5p for S, we see 

, . b 2 - a? . x - 1 / 6 - o 

1*2)12 g 2^ ^rn" I 4 iH-mlcrs + -Ttm2<Ti 

m odd 

El (b — a 



W 
in even 



and 



(^1)22 ~ - X] 7777^™^ + ^ 

m od 

From (14.21) it then follows that 



mil 

Ti" h " ' J Tt"' 

m odd m even 



m=0 



where /3 = ^ and 



ft = 2R llas - J—iR 12ai + 2(6 + a) ^ llgl . (4.9) 
— a — a 

Our final task is to further evaluate the right-hand side of (|4.9p . As in 
[14j . we have that Ai is meromorphic in a neighborhood of b with a pole in 
6. Indeed, if we write 



(5{zY 



(z-b)-^B m (z-br, B = —— (4.10) 



m " ' to ' ' U 2irh(b) ' 

and use (|3.22p . then we find for z in a neighborhood of 6, 

5-Bi 75 , . A 1 



A l( Z ) = ( -777 (°3 - tO-l) + t,,/, \ (^3 + M 



144 v ' 144(6 -a) v ' ) z-b 



144 v {z-b) 
Similarly, for z in a neighborhood of a, we have 

/3(z) 2 



( z -a)- 2 Y, A m(z-a) m , A = —— (4.12) 
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and 



A lO) = ( ~7T7 (°3 + i(T\) - , AA ]t° , (0-3 - iffi) 



144 v ' 144(5 -a) v 7 z-a 



As in [14] we have that R\ (z) for z G C\U U J7 is equal to the sum of the 
Laurent parts of (|4.11j) and (|4.13|) . Expanding Ri(z) as 2: — > 00, we then get 

= R u - + + o(\) as z ^ 00, 



where 



z z 



^ n = " 144 (CT3 + Wl) " 144(6 - a) ^ ~ 
5-Bi 7S 
- 144 (a3 ~ mi)+ 144(6-a) ((73 + ^ 
5aAi 7ayl , . s 

12 = ~ HT (as + mi) " 144(6 -a) (CJ3 " iai) 
5bB, 7bB . 

-T4T (a3 - mi)+ 144(^) (fT3 + mi) 

" 144 fo + itfi)- — (as-mO. 



Thus 



_ 5(Ai + ^i) 7(^o - go) 
^ Wa 144(6 - a) ' (414) 

^ = "* 144 + 1 144(6 - a) ' (415) 

.5(aAi-6Bi) 7(ai4o + 6B ) .5(4, --Bo) 
^ 12 - = 144 + 1 144(6 -a) ~ * 144 ' (416) 

Inserting (j4.14fl — (j4.16H into (|4.9f) . we find after straightforward calculations 
that A\ and -Bi fully disappear and that (|4.9p reduces to 

_ £ - A) 
Pl 3(6 - a) ' 

Using the explicit formulas for Aq and Bq given in (|4,10p and (|4.12p . we 
arrive at ()1.5|) , which completes the proof of Theorem 11.11 □ 
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